DISPERSIVE ESTIMATE FOR THE SCHRODINGER EQUATION WITH 

POINT INTERACTIONS 



PIERO D'ANCONA, VITTORIA PIERFELICE, AND ALESSANDRO TETA 

Abstract. We consider the Schrodinger operator in R 3 with N point interactions placed at 

Y = (j/i, , Uff), yj £ R 3 , of strength a = (ax, • • • >cin)> a j € K. Exploiting the spectral 

theorem and the rather explicit expression for the resolvent we prove the (weighted) dispersive 
estimate for the corresponding Schrodinger flow. 

In the special case N = 1 the proof is directly obtained from the unitary group which is known 
in closed form. 



1. Introduction 

In this paper we study the time dependent Schrodinger equation in R 3 with a finite number of 
point interactions and, in particular, we shall prove a dispersive estimate for the solution. 

At a formal level the Schrodinger operator with point interactions can be written as 

N 

(1.1) ff = -A + £>A, 

i=i 

where S Vj is the Dirac measure placed at yj £ R 3 and the parameters /ij are coupling constants. 
As a matter of fact the Dirac measure in R 3 is not a small perturbation of the Laplacian, even in 
the sense of quadratic forms. 

As a consequence a self-adjoint Hamiltonian in L 2 (R 3 ) cannot be defined as the sum of a kinetic 
plus an interaction part and this explains why (1.1) is only a formal expression. 

In order to obtain a rigorous counterpart of (1.1) one considers the following restriction of the free 
Laplacian 

(1.2) H=-A, D(H) = C^°(R 3 \Y) 

where Y = (t/i, . . . , Vn)- The operator (1.2) is symmetric but not self-adjoint in L 2 (R 3 ) and, 
obviously, one possible self-adjoint extension is trivial, i.e. it coincides with the free Laplacian 
H Q = -A, D(H ) =H 2 (R 3 ). 

Using the theory of self-adjoint extensions of symmetric operators, developed by von Neumann 
and Krein, one can show that the operator (1.2) has N 2 (non trivial) self-adjoint extensions which, 
by definition, are all the possible Schrodinger operators with point interactions at Y (for a com- 
prehensive treatment we refer to the monograph [1]). 

Any such extension can be considered as a Laplace operator with a singular boundary condition 
satisfied at each point yj G Y. 
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In the following wc shall only consider the case of local boundary conditions which are more relevant 
from the physical point of view. More precisely, we shall restrict to the self-adjoint extensions H a y 
parametrized by a = (a\, . . . , ajv), aj € M, and corresponding to the singular boundary condition 
at Y 



(1.3) 



lim 



dr.i 



0. 



3 = 1, 



In the next section we shall give the precise definition and the main properties of the Hamiltonian 

We only notice that the domain of H a ^y contains functions with a singularity of the type \x — yj | _1 
at each point yj e Y and this explains why in the boundary condition (1.3) the behaviour of the 
function u near yj must be regularized. 

We also recall that the physical meaning of the parameters ctj is connected with the scattering 
length of the scattcrer in yj and their relation with the formal coupling constants fXj introduced 
in (1.1) can be understood via a suitable renormalization procedure (see e.g. [1]). 

In the following we shall be concerned with the unitary evolution group generated by H a y 



(1.4) 



JtH a 



which gives the solution of the Cauchy problem 



(1.5) 



iu t + H a ^ Y u = 0, u(0,x) = f(x). 



It is well known that for the free Schrodinger group in dimension three the following dispersive 
estimate holds 



(1.6) 



\e- itA f\ 



< 



C 



L oo( R 3) — £ 3 /2 



The estimate (1.6) can be generalized to smooth perturbation of the Laplacian, i.e. the l.h.s. can 
be replaced by -P a ce J *^ A+y - 1 /! where V is a smooth real-valued potential satisfying some decay 
condition at infinity P ac is the projection onto the absolutely continuous spectrum of —A + V and, 
moreover, it is assumed that zero is not an eigenvalue nor a resonance for —A + V (see e.g. [3], [4] 
[6], [7] and the references therein). 

Here we want to extend the validity of the dispersive estimate to the case of the Schrodinger 
equation with point interactions. 

Due to the presence of the (unavoidable) singularity at the points where the interaction is placed, 
we are forced to introduce the following weight function 



N 



(1.7) 



w 



Moreover, for z e C, wc define the matrix 



(1.8) 



\Pa,Y{z)]j,t 



a 3-^j s i,e - Gz(yj-ye) 



N 



3,t=l 
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with 

{ e iz\x\ 
x = 0. 

The role of the matrix r aj y(z) for the characterization of the main properties of H a y will be 
explained in the next section. Here we only notice that our basic assumption is the invertibility 
of T a ,y(n), for yu £ [0, +oo). Such assumption in particular implies that zero is not an eigenvalue 
nor a resonance for H a y- 

In the special case N = 1, i.e. for the Schrodinger operator H atV with a point interaction in y £ R 3 
of strength a £ K., the assumption simply means q / 0, which is precisely the condition for the 
absence of a zero-energy resonance. 

Our main result is the following theorem. 

Theorem 1.1. Assume that the matrix r Qi ^(/z) is invertible for ji £ [0, +oo) with a locally bounded 
inverse. Then the following dispersive estimate holds 

(i.io) KV'^.-p ac /|| ioo(R3) < ^ IK /Hlw 

for any f £ L 2 (R 3 ) such that w ■ f £ L^M 3 ). 

In the special case N = 1, estimate (1.10) holds for all a ^ 0; moreover, when a > the projection 
P ac can be replaced by the identity. Finally, in the resonant case a = we have the slower decay 
estimate 

(!- n ) |KV«^/|| ioo(R3) < £s\\w • /|| Ll(R s) 

We remark that the result for N — 1 is more detailed due to the fact that the unitary propagator 
is explicitely known (see e.g. [5], [2]) and then the estimate can be obtained by a straightforward 
computation. 

The rest of the paper is organized as follows. 

In section 2 we recall the precise definition and the main properties of the Hamiltonian H a y. 

In section 3 we give the proof of the theorem in the special case N = 1, exploiting the knowledge 
of the kernel of the Schrodinger group. 

In section 4 we give the proof in the general case using the spectral calculus and the explicit 
expression for the resolvent of H a ,Y- 

2. The Schrodinger operator with point interactions 



In this section we review the definition and some basic properties of the Schrodinger operator 
with point interactions, referring to [1] for more details. We start considering the perturbation of 
—A by a zero-range (singular) potential supported at the origin of IR 3 . As we pointed out in the 
introduction, such perturbation is by defmiton a non trivial self-adjoint extension of the symmetric 
operator (1.2) when N = 1. 

We recall that the trivial extension —A is the (self-adjoint) Laplace operator in L 2 (M 3 ) (with 
Lebesgue measure), with domain D(-A) = H 2 (R 3 ) = T^ 2 - 2 (IR 3 ), where W p ' k (M. 3 ) denotes the 
Sobolev space of all functions belonging to L P (K. 3 ) whose weak derivatives of order smaller or 
equal to k belong to i p (M 3 ). 

Using the theory of self-adjoint extensions, one can show that the non trivial self-adjoint extensions 
H a , parametrized by a £ R, can be completely characterized. 



(2.1) Mx)=Mx) +[<*-¥-) MO) 
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Indeed, the domain V{H a ) consists of all elements tp <G i 2 (M 3 ) of the type 

:\x\ 

An J ^ '~ K ' ^n\x\ 

where <p z <G H 2 (R 3 ) and 3z > 0. The decomposition (2.1) is unique and with ip e T>(H a ) one has 

(2.2) {H a -z 2 )i, = (-A-z 2 )<j> z . 

We notice that the trivial extension is recovered in the limit a — > oo. 

A remarkable property of the Hamiltonian H a is that the integral kernel of the resolvent R a (z) = 
(H a — z)~ x , 3z > 0, can be explicitely computed by the krein's formula. In fact 

(2.3) (R a (z 2 )f)(x) = (R (z 2 )f)(x) + (a- l -f] ' [ ^\f(y)dy. 



4irJ 4ir\x\ J m 4ir\y 



where the free resolvent is given by 

Jz\x-y\ 



(R (z 2 )f)(x) = / ±- -f(y)dy, 

Jr3 A-k\x - y\ 

By an elementary computation we obtain that for any Ael and e > 

f p ±iVX^\x-y\ 

(2.4) (R (X±is)f)(x)= — e-^-y\l 2 ^f{y)dy, 

J R 3 \-k\x - y\ 

where 

(2.5) A , = A±(j!±f!)^ >0 . 

From the last inequality it is easy to derive the limiting absorption principle for the free resolvent, 
i.e. for A > 

p iyf\\x — y I 

(2.6) (R (X + iO)f)(x) = 6 — -f(y)dy. 

J R3 A-k\x - y\ 

We can also deduce the limiting absorption principle for the resolvent (2.3) 

(j\/a\ 1 p-'v'Alxl r p i\/~\\y\ 
a -in-) ~lnwL^\yT f{y)dy - 

The spectral properties of H a are easily derived from (2.3). The continuous spectrum of H a is 
purely absolutely continuous and covers the nonnegative real axis [0,oo) while the point spectrum 
is empty if a > and {(— 47ra) 2 } if a < 0. The normalized eigenfunction associated with the only 
negative eigenvalue is given by 

e A-na\x\ 

(2.8) ^Jx) = V-2a , , , a < 0. 

M 

For a = the Hamiltonian has a zero-energy resonance. 

In analogous way we introduce the properties of the Schrodinger operator H a y with point inter- 
actions located at Y = (j/i, . . . , y^) with strength a = (ai, . . . , a^). 

The domain 2?(_ff Qi y) is the set of all functions ip G i 2 (R 3 ) of following type 



(2.9) ^) = ^WlVa ir - - x£B 3 \F, 
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where 

JV 

a j = 5^[Ta,y(«)]^Vz(w). j = l,.-.,N, 
1=1 



<f> z € H 2 (R 3 ), 3z > and the matrix T a ^y{z) has been defined in (1.8). The decomposition (2.9) 
is unique and with ip <G V(H a ^y) we obtain 

(2.10) (H aiY -z 2 W=(-A-z 2 )<j> z . 

It is an easy computation to verify that each element of V{H a ^y) satisfies the local boundary 
condition (1.3). 

Also in the TV-centers case the resolvent of these operators is explicitly given by Krein's formula. 
In fact R a ,y(z) = (H a ,Y — z)^ 1 is given by the following formula for 3z > 

(2.11) (R a .y{z 2 )f)(x) = (R a (z 2 )f)(x) + J2 Fa, y (^)] J -l f "' a 7'| / *T~ V '\ mdy 

fj± x ^A x - Vj\ Jm ^\y - vi\ 

Moreover we deduce the validity of the limiting absorption principle, i.e. for A > 
(2 - 12) 

J-L _ e -iV^\x-yj\ f „iV*\y-ye\ 

(R a . Y (X + i0)f)(x) = (R (X + i0)f)(x) + X: ]T atY (iy/X)]j} ^-—^ f(y)dy 

We conclude summarizing the spectral properties of H a _y. The essential spectrum of H a is purely 
absolutely continuous and coincides with the real axis [0, oo), indeed the singular continuous spec- 
trum is empty. Moreover, the operator H a y has the point spectrum included in (0, — oo), i.e. there 
are no positive embedded eigenvalues; in particular H a ^y has at most N (negative) eigenvalues 
counting multiplicity. In addition, the eigenvalues can be determined computing the zeros of the 
determinant of an iV x N matrix, i.e. 

z 2 S cr p (H a , Y ) detr Q ,y(z) = 0, 

and the multiplicity of eigenvalue z 1 equals the multiplicity of the eigenvalue zero of the matrix 

T a .y(z). 

3. Proof of the result for N = 1 

In the case of a single point interaction considered in in the second part of Theorem 1.1, the proof 
is quite easy. We believe it is of some interest, both because the proof is much easier, and it is 
possible to give a complete description of the dispersive behaviour of the propagator, including the 
case when a resonance at zero occurs. The main advantage is that, as proved in [5] (see also [2]), 
it is possible to give an explicit representation of the Schrodinger propagator. If we denote by 

\x\ 2 /4it 

(3.1) S(x; t) = 6 , t>0, x,yeR 3 

the free Schrodinger propagator in R 3 , i.e. 

(e- itA f)(x)=S(x;t) * x f(x), 
the kernel of the propagator of e ltHa can be written as follows: 

1) for a > 

1 f°° 

(3.2) S a (x, y; t) = S(x -y;t) + ---- / e-^ as (s + \x\ + \y\)S(s + \x\ + \y\;t)ds 

\ x \\y\ Jo 

2) for a = 

lit 

(3.3) S (x, y; t) = S(x -y;t) + t^S(M + \y\;t) 

\ x \\y\ 
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3) for a < 
(3-4) 

S a (x, y, t) = S(x -y;t) + ^ a (x)^ a (y)e lt ^ 2 + ^ [°° e^ as (s - \x\ - \y\)S(s - \x\ - \y\;t)ds 



mm Jo 

Then we have the representation 



(e itH «f)(x)= f S a (x,y;t)f(y)dy. 

JR3 



Notice that in the case a < we can easily distinguish the continuous part from the standing wave: 
indeed, the continuous part can be expressed as 

(3.5) (P ac e itH <> f){x) = (e- UA f)(x) + f f°° -^-e 4 ™ s ( S - \x\ - \y\)S(s - \x\ - \y\;t)dsdy. 

Jr. 3 Jo \ x \\y\ 

The standing wave is the function 

Ci^aOe"* 4 ™) 2 , Ci=/ * a {v)f(v)dy 

Jr 3 

which can be written explicitly using (2.8) as 

(3.6) (-2a)-i-r- / -rT-f{y)dy. 

fi Jr 3 \y\ 

Consider the case a > first. From the trivial estimate 

(3.7) |5(« + |a:| + |y|;t)|< 



Consider the case a > first. From the trivial estimate 

C 

t 3/2 

and 



(3.8) 



/ e - 4 ™ s (s + \x\ + \y\)ds<C(l + \x\ + \y\) 
Jo 



and using also the standard dispersive estimate for the free solution, we obtain immediately the 
estimate 

S a (x,y;t)f(y)dy < tSj / fl + 1 ~ ) \f(v)\dy. 

1 1 Jr 3 V \ x \\y\ / 

Since we have 

1^1 12/1 V fi/ V 

we easily conclude that 

\( w -'e^f)(x)\ < ^ J^w( y )\f(y)\dy 

as claimed. 

In the case a < 0, we can estimate the continuous part (3.5) exactly in the same way, obtaining 

\{w^P ac j tH "f){x)\ < J^w(y)\f(y)\dy. 

Finally, when a = we can see the effect of the resonance at zero in the slower rate of decay: from 
(3.3), proceeding exactly as before, we obtain 

l^-V"^ /)(*)! < -gj / w(y)\f(y)\dy 
t L/z Jrs 

since an additional factor t is present in the propagator. 
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4. Proof of the result for N > 1 

We shall now prove the main part of Theorem 1.1 concerning the general case N > 1, where we 
do not have an explicit formula for the propagator. Thus we resort to the spectral calculus and we 
represent the (continuous part of the) solution as follows: 

r+oo 

(4.1) P ac e* H <"*f = / e itx SR a , Y (\ + iO)fd\. 

Jo 

Recall now the integral expression (2.12). The first term reproduces the free solution e~ ttA f for 
which we already know the decay estimate. Thus we are reduced to estimate the integrals 



-dydX 



where 

(4.3) Cjtin) = {Air)- 2 [T a , Y (ii)]-j} for M > 0, j,£ = l,...,N. 

To make the following computation rigorous, we introduce a cutoff function. Let i/)(s) be a non- 
negative function in C^°(R) equal to 1 on [0,1] and vanishing on [2, oo); it is not restrictive to 
assume that tp is an even function ip(—s) — ip(s). Then we can approximate / with 

(4.4) I M = [ + ° ° f e^sfe-'^'-^j^-ricjeiVX))-, ^ ^{V\/M)dyd\ 

Jo Jrs v ' \x-yj\\y-yi\ 

asM^ +oo. After the change of variables A = fj, 2 this can be written also as 

(4.5) Im=[ +C ° I 2^^(e- i ^ x -^e i ^-^c j ^))- -Mn/M)dydn 

Jo Jm v 1 \x-yj\\y-yi\ 

If we can prove a dispersive estimate for Im which is uniform in M, this will also give an estimate 
for /. Since cje(fi) is bounded near \i = by assumption and c^(0) is real, while 

cj( e -vk-wl e vlw-w«l) = S m( M (|a; - y e \ - \y - yj\)) 

vanishes for \i = 0, we can integrate by parts with respect to \i and we obtain 

(4.6) I M = \ [ +C ° I e^'dpS (e-^-^e^-^cM^/M)) rdyd^i. 

t Jo jk 3 v J \x- yj\\y - y e \ 

For reasonable data / the two integrals can be swapped using Fubini's theorem and we arrive at 

(4.7) I M = 7 / -. ^ r / °° e*"' 8^ (e-^-^e^-^Cje^iii/M)) dfxdy. 

tJ R s \x-yj\\y-ye\ 7 v ' 

The core of our proof will be a suitable estimate of the inner integral in the variable \x\ expanding 
the derivative we obtain three terms 

" e^'d^ (e-^-^e^-^cM^/M)) rf/i = h + h + h 



f 

Jo 



10 

and precisely 

f + OO 



r+oo 

(4.8) h = (\y - ye\ -\x- yj\) J e'*" 3? (e^-^-^-^c^M/V^)) d/x, 

(4.9) I 2 = J Q^e-^-^e^y-^c'^^/M^dfj,, 
and 

(4.10) h = e^ 2 3 (^-^-^e^y-y^cMii'iii/M)^ dfi. 
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Consider the first integral I\ , or rather the integral 

(4.11) h = J e Utl 0? ye i ^ y ~ y ^~ i ^ x ~ y ^Cjt(fj,)'ip(iJ 1 /M)J dfi. 

Inspired by a clever idea of Rodnianski and Schlag ( [4] ) , we remark that an integral of the form 

/■+OO 

u(t, A)= e'*" 3? (e iA »c je (iJ,)ip(ti/M)) dn 
Jo 

is a solution of the one dimensional Schrodinger equation idtu — d\u. Thus we can apply the 
classical dispersive estimate 

\u{t,A)\ < Cr 1 ' 2 J \u(0,A)\dA 

which in our case gives 

C 



h 



- ■ J |Ji(A)|dA 



Vt 



(4.12) 
with 

r+oo 

(4.13) Ji(A)= / K(e iA i*Cjt(jJ.)il>(ji/M))dii. 

Jo 

The same argument, used for I 3 , shows that 

(4.14) |/ 3 | < • J \J 3 (A)\dA 
where 

(4.15) MA) - J o + °° 9 [e^ A c lt (^'^/M)^j d/x. 

A similar argument can be applied to I^. First of all we recall that the matrix r Q .y(/i) for 
[i G [0, +oo) has the form (see (1.8)) 

(4.16) r a>y (/i) = -^-1 + diag[ai, ...,a N ] + A(ji) 
where A(fi) has coefficients 

e iv\yj-yt\ 

(4.17) ^(M)k=0, [A(fi)] je = - - iovj^t. 
Hence the derivative 9 /j r Qj y(/x) is simply 

(4.18) [d»T a , Y {i i )] je = -±eWy>-y<l 

Since 

d^a^Yipr 1 } = r a ,Y(py 1 (d fJ ,T a . Y (iJ.))T a .Y(py 1 , Cjl{n) = [r^w^r 1 ].^, 

we obtain the representation 

N 

(4.19) C ^ (/U ) = _J_ cjkMe^-^CkM 

k% ,k 2 = l 

If we plug this into Ii we obtain 

(4.20) J 2 = — V / e^ 2 sfc jfe J/i)c te< (M)e ip| ^ w| -^ |x -w |+ ^ lw i-^^(/i/M)) d/z. 

4 vir=i 7 ° v 7 

The same argument used above for I\ and I 3 gives that 

(4.21) |7 2 | < ^ • / I J 2 (A)|dA 
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where 

r+oo 

(4.22) MA)= Bfe^CjfeWcwWfWM))^ 

Jo 

We shall now prove the estimate 

(4.23) J (IMA)] + \J 2 (A)\ + \J 3 (A)\)dA < C 

for some constant independent of M, from which the Theorem will follow immediately. To this 
end we shall need to study the coefficients Cjt(fi) closer. 

Lemma 4.1. Let T a ,Y(p) be the matrix defined in (4.16), (4.17) and assume it is invertible with 
inverse locally bounded for fi G [0, +oo). 

Then the coefficients Cje(n) of the inverse matrix r ai y(/z) _1 satisfy the following properties: 

(i) The coefficients are holomorphic on a neighbourhood of the positive real axis; their derivatives 
satisfy the estimates 

(4.24) \cM\ < cw 1 , \ c ' m\ <c(^)-\ \c';m\<c^)-\ M >0 

where we used the notation (/i) = (1 + /i 2 ) 1 / 2 . 

(ii) The coefficients can be written as 

(4.25) Cji(^) = ^iriSji^y 1 + dje(n), 
or equivalently 

(4.26) Cji(p) = 4:TTiSjHJ,(iJ,y 2 + dji(fi), 
where the holomorphic functions djg,dj t g satisfy the inequalities 

(4.27) \d ol \ + \d' ]l \ + \d'; i \<C^)- 2 , \d ]l \ + \d' ]l \ + \d'; i \<C^)-\ M >0. 

Proof. It is clear from the definition that Cji(n) are holomorphic (as soon as they are defined); 
moreover, by formula (4.19) proved above we see that the second and third estimates in (4.24) are 
immediate consequences of the first one in (4.24). 

In order to prove the first inequality in (4.24), it is sufficient to recall formulas (4.16), (4.17), i.e., 

r Q ,r(M) = C fil + D + A(fi), 

with Co = — i/An, D a constant diagonal matrix, and A(fi) a matrix with bounded coefficients. 
From this expression it is clear that the entries of the inverse matrix r ai y(/i) _1 are of the form 

±(C oM ) w - 1 +P( M ) 
CjeW ~ ±(Con) N + Q(») 
where P, Q are functions of \x of order N — 2 and N — 1 respectively: 

\P(JJ)\ < C([i) N ~\ \Q(jj)\ < CM*- 1 . 

Taking into account the assumption that cje(fi) are locally bounded functions, estimate (4.24) 
follows easily. 

In order to prove the asymptotic expansion (4.25), we write 

I = T^Yifi^T^Yifi) = T^yi^iC^I + D + A(jj,)) = Co/J^yOu) -1 + T^y^Y^D + A(ji)] 
which implies (Cq 1 = Airi) 

r^viv)- 1 = Airi^L - [D + A(n)] ■ r^Y^yHnin- 1 
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(it is sufficient to prove the estimates for \x > 1 since we already know that the functions are 
smooth near 0). Notice that the last term is bounded by C(/i)~ 2 . Now, using the elementary 
identity 

1 _ 1 1 

M ~ (1 + M 2 ) 1/2 + Ml + A* 2 ) 1/2 [^ + (1 + M 2 ) 1/2 ] 
and estimate (4.24) already proved, we easily obtain (4.25) and the estimate 

\d je (u)\ < C( M )- 2 . 

The remaining estimates on the derivatives of djt follow immediately from (4.25) and (4.24). 
The second expansion (4.26) is proved in an identical way; indeed, 

1 n 1 

M ~ 1 + A* 2 + M(l + M 2 )' 

□ 

Remark 4.1. The reason for the two different asymptotic expansions in the Lemma is the following: 
the main term decays like /x _1 and must be treated with care, using explicit formulas for the Fourier 
transform. Since the Fourier transform here appears as a sinus or cosinus transform only, it is quite 
useful to have its expression both as an even and an odd function of \i. The remainder terms are 
harmless since they give rise to integrals which can be easily estimated uniformly in M. 

We are ready to estimate the L 1 norms on R of the three quantities Ji(A), J2(A), J%(A). Consider 
the quantity Ji(A); using the expansion (4.26) we have 

(4.28) Ji(A) = J" (A) + J\(A) 

where 

p+ca r+co 

J?(A) = 4tt / 3? (e lAf "4TTi^)- 2 ^/M)) d/x = -(4tt) 2 / sm(Afj,) l j,{fj,)- 2 tp( t j,/M)d f j, 
Jo Jo 

and 

r+ca 

J b M) = J 3? (e^^W^WM)) d/x. 

The remainder term J b (A) can be estimated directly. Indeed, by (4.27) it is easy to see that J b (A) 
is bounded 

r+oo 

\Ji(A)\ < C / (n)- 2 dii<C 
Jo 

with C independent of M,A. Moreover, writing e tA ^ = (iA)~ 1 d l _ l e lA ^ we have 



1 f +oa / \ 



but Cji(Q) is real, so integrating by parts we get 

1 r +cc / \ 1 f + °° / 1 \ 

4(A) = -jj o 9 (e^d'^^/M)) d^--J^ 3 [e^djfaWWM)—} d/x. 

As before, using (4.27) we see that both integrals are bounded by a constant C independent of 
M, A and this implies 

\AJl{A)\ < C. 

Finally, we write e tA ^ = (iA)~ 1 d lx e lAl1 one more time in the last formula: 

A{A) = -1 ^ + °° 3? (d^di^^/M)) d M + £°° n (a.e^^Mf'WM)^ d/x. 

We integrate again by parts; this time we obtain a boundary term 
(4.29) ±-m' je (0)m + -^ZtdjdOWWjj 
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and a few integrals 



(4.30) — jf S^^WWM))^-! R^dJ^/i^WM)— Jd/i 

(4.31) +-1 ^ + °° 5ft ^"^(^jfWMjlj d/x. 



Applying as above the estimates (4.27) of the Lemma we obtain 

\A 2 j\(A)\ < C" 
for a constant independent of M, A. Then we conclude 



J \j\{A)\dA < J C(l + A 2 )- x dA < C 



for some bound C uniform in M. 

Let us now consider the main term Jf (A) in (4.28). It is a standard sinus-transform, which can 
be rewritten as a Fourier transform as follows (recall tp is an even function) 

/+oo 
e ^7T^W^)^ = C ■ {T{ 9l )*T(g 2 )) (A) 
-oo 1 + M 

where we have used the convolution theorem and we have denoted 



(4-32) 9l {p) = ff2 (/i) - VWM) 

1 + M 

Since 

(%))(4)=cA e -|A| 

and 

(^(32)) (A) = M(f^)(M) 

we obtain 

KIlw < c||H9i)IIl 1(R ) IIH92)IU 1(R ) < c 

for some constant independent of M. In conclusion we have proved the bound || JiIIl^r) < C, 
uniformly in M. 

The corresponding estimate for Js(A), see (4.15), is almost identical; the only modification is the 
use of the explicit formula for the Fourier transform of (fJ,)~ , which is 

H{H)- 1 ) = CK Q {\A\). 

Here K (s) is the modified Bessel function of order 0, whose behaviour is the following: 

\K (s)\ < qiogs|, 0<s<l; \K (s)\<Cs-h- s , s > 1. 

In particular, it is clear that J 7 ^)^ 1 ) belongs to L 1 (M) which is what is needed in our computation. 
On the other hand the estimate for J2(A) (see (4.22)) is immediate since the integrand decays as 
(^i)~ 2 . This concludes the proof of (4.23). 

Recalling now (4.8), (4.9), (4.10) and (4.12), (4.21), (4.14) we have 

\h\ < ^(\V-Vl\ + \X-Vi\), \M + \h\ < ^ 

with C independent of M and of x. Thus, recalling (4.6) and taking the limit as M — > oo, we 
obtain the following estimate for the integral I in (4.2): 



\i\<£rJ \f(y)\ 1+ ! y - y i + lx - yjl dy 
t 3/2 Jr3 \y-ye\\x-yj\ 
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Since Pa C e ttHa Y f is the sum of the integrals / for all couples j, £ = 1, . . . , N, plus the free wave 
e ~ ltA f, w e obtain 

The inner sum gives 



where w(a;) is the weight function (1.7). 
In conclusion we have proved 

wixr^P^*^ f){x)\ <-^J \f(y)\w(y)dy 

as claimed. 
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